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Abstract 

Vn  analytical  imestigation  of  ilu*  ihiumu  performance  ol  .1 
centrifugal  lilt  Ian  was  conducted  to  provide  an  explanation  tor 
the  behavior  which  occurred  when  tlu  Ian  was  Mihicctcd  to  a 
vaivinji  backpressure  Analysis  ot  test  data  shows  that  the  Ian 
response  M’  can  be  represented  bv  a lust-ordet  lag  svstem 
An  in  depth  analvsis  ol  the  various  elements  ot  the  lilt  tan 
system  revealed  that  the  inertia  of  the  air  within  the  fan  was  the 
primary  contributor  to  the  observed  Ian  behavior  I he  analysis 
further  showed  that  variations  in  Ian  speed  would  not  pioduce 
the  behavior  measured  in  the  test  and  that  the  response  due  to 
the  compliant  pioperties  of  an  within  the  tan  occurs  at  a 
frequency  that  is  much  higher  than  the  frequency  i.uige  ol  the 
test  v time-domain  digital  computer  program  has  been 
developed  which  integrates  the  rate  of  change  of  Ian  floyv  with  a 
varying  backpressure  Good  correlation  is  exhibited  between 
test  data  and  the  computet  predictions  at  all  frequencies 

Nomenclature 

A Area  ol  flow  channel,  ft* 

A,  Area  of  i*"  section  ol  plenum  wall,  ft* 

U Compliance  of  air  in  fan.  It'  lb 

C|>  Discharge  coefficient 

I ; Inverse  ol  fan  characteristic  slope,  it'  s/lb  If 

f I lequencv . II; 

I Inertia  ol  rotating  system.  11  Ibf  s* 

1 (-i*is 

l I net  lance  of  air  in  fan.  Ib-s*  It' 

m Mass  of  an  in  the  fan.  Ibm-s*  It 

N Shaft  rofafional  speed.  rpm 

P Pressure,  lb  tt: 

0 I an  flow  tale,  ft  * s 

O,  Uncontrolled  plenum  floyv  due  to  leakage,  ft' s 

Op  Controlled  How  from  plenum,  ft' s 
S Distance  along  flow  path,  ft 

1 Torque.  Ibf  It 

t Time,  s 

\p  Plenum  volume,  ft' 

\ Velocity.  It  s 

/,  IVflection  of  i**'  section  of  plenum  wall,  ft 

fi  Mass  density  . Ibm-s*  ft4 

u.’  Rotational  velocity,  rail  s 

Subscripts 

\IM  Ambient  conditions 

c Corner 

1 1 .in  ' 

OnpywfigM  Akwnw  Inalitut*  of  AhinmuIiti  «nri 
AifronculKS.  Inc.  ItTl.  All  fight*  react  veil 
N*C«ni1iic«rt  with  pa*  mwtnn 


I Inlet 

I I cakagc 

M Motor 

P Plenum 

S Static 

I Total 

Introduction 

In  Surface  I Meet  Ships  r SI  Si  and  \u  Cushion  V ehicles 
i ACA  ».  the  lift  fan  characteristic  pressure  versus  flow  late  curve 
is  known  to  have  a significant  effect  on  the  vertical  motions  ot 
the  platform  I xposurc  of  the  ciew  to  adverse  motions  results 
in  degraded  human  performance  and  teduced  overall  system 
ettectiy eness  1 Analysis  techniques  to  predict  the  motions  of  SI  s 
and  A(’\  s span  the  spectrum  from  simple,  lineai.  single  degree 
ol  tieedom  models  to  elaborate,  time-domain,  multi-degree  ot 
tteedom  models  In  all  cases,  the  primary  xctttcal  force  arises 
from  the  cushion  pressure  acting  on  the  wel  deck  The  cushion 
piessure  is  derned  m these  analytical  models  In  integrating  the 
time  rate  ot  change  ot  the  mass  of  air  in  the  cushion  plenum 
which  in  turn  is  determined  using  the  conservation  ot  mass 
equation  The  late  at  which  an  is  dehxered  from  the  fans  to  the 
plenum  is  usually  determined  using  the  steady  state  fan  chatac 
tenstic  in  the  form  0 ft  pi  This  assumes  that  the  tan  behaxior 
is  the  same  under  steady  state  and  dynamic  conditions 

In  ordei  to  inxcstigatc  the  behaxior  of  hit  tans  in  unsteady 
flow  conditions,  a test  rig  was  constiuctcd  at  the  D.txidW  l ax  lor 
Naxal  Ship  Reseat ch  and  Dexelopment  l entei  iDINSRDO  I his 
lilt  fan  exaluation  ng  features  an  S'O  ft  plenum  and  an  exhaust 
duct  containing  a flow  control  xalxe  d ig  1' 


I ig  I - Artists  Sketch  of  the  Tan  I valuation  Rig 
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Ilk*  Ian  delivers  pressurized  an  i nit*  the  I'K- mu  in  \ppropnaic 
instrumentation  is  hualed  it  the  fan  inlet  aiul  dis*haige  to 
measure  the  Mow  iate  and  |>(fttuiv  iiw  across  the  tan  \ii 
evaluation  was  *onductc*l  using  a '4  in  diameter  m**del  ol  the 
si  n IOOH  cushion  tan  representative  ot  high  eMuietuv  back 
waul  swept  an  toil  till  li  V » *entnfugul  tans  Data  horn  this 
evaluation  showed  that  uiulet  dvnatnu  conditions  the  How  iate 
change  lagged  the  ptovsure  change  in  su*h  a mannei  that  the 
pieviute  versus  Mow  iate  *haia*  tetisti*  iu*  longer  followed  the 
sleadv  state  curve  hut  instead  when  plotted  on  the  standaid 
pieusutc  versus  Mow  iate  plot  exhibited  loops  l ig  piescnts 
examples  ot  the  piessuie  versus  flow  iate  hvps  that  vveie 
obtained  during  the  evaluation  ot  the  sl  s 100H  scale  nunlel 
cushion  tan  these  loops  are  superimpose*!  on  the  steals  state 
fan  *haia*  tenstl*  t he  data  show  the  variation  of  the  loop  shape 
with  the  !rcs|uen*v  of  the  back  pievsurv  change  Ihev  alsv* 
illustiate  the  steepening  ol  the  die* five  sU*pe  of  the  fan  vhaia* 
tenstu  uletemuneil  hv  the  orientation  of  the  1114101  avis  of  the 
“ellipse"  tonne*!  hv  the  loop!  as  the  H*s|uen*v  is  increase*!  \ 
lift  sv  stein  designer  could  tight  lull)  ask  what  tan  chatactenstu 
shoulvl  he  select*'*!  lot  tvpical  dvnatnu  c**n*!iti**ns 


» ig  : SI  S IOOH  Seale  Model  l lit  t an 
Ih na line  lest  Results 


I his  paper  presents  the  results  **t  analvses  comlucfo*!  on  the 
data  obtained  in  the  evaluation  ot  the  scale  m*\!el  centrifugal 
lift  fan,  and  a method  for  predicting  fan  behavior  m a dvnami* 
environment  V computer  simulation  based  on  the  method  is 
verified  using  experimental  data  from  the  tests  of  the  Si  S IOOH. 
model  fan,  as  well  as  a I 4 scale  model  of  the  'kSl  S lift  Ian 
I he  analytical  nunlel  of  tan  dynamic  performance  should  In- 
valuable in  simulations  of  vertical  motion  of  an  SI  S oi  U \ in 
a seawav 


Analysis  ot  Fan  Dynamic  Data 

l ime  histones  of  total  and  static  pressure  rise  and  an  How 
at  frevtuencies  of  0 s H/  and  (*  11/  are  shown  in  I ig  * These 
data  are  illustrative  of  data  obtained  during  the  test  ot  the  '4  in 
diameter  S|  S liH>R  nunlel  cushion  fan  Data  were  collected  at 
fan  speeds  ot  _'<c*0  and  ' rpm  NUn-ie  the  wave  forms  were 
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sultuu-ntlv  sinusoidal,  tusjueiuv  lesponse  mtoimation  in  teims 
**t  the  amplitude  ratio  and  phase  angle  between  tan  Mow  aiul 
total  pievsuie  use  wete  evaluated  I xammation  ol  the  tnnjueiuv 
resp**nse  intonnaii**n  indicates  that  the  amplitude  lain*  decreased 
and  the  phase  angle  me  teased  as  the  excitation  inxjuencv  was 
in*  leased 

I he  response  toi  a tuspieiuv  o!  0 s and  I 0 II  exhibited 
very  large  variatu*ns  m M**vv  and  piessuie.  induating  extremely 
nonlinear  opeiation  I he  operating  p**mt  vaned  fiom  neatly 
untestneted  Mow  t**  a Kick  Mow  eoiulition  \t  the  high  Mow 
pot  lions  ot  each  cycle  the  Mansion!  lesultmg  tn*in  the 
nonlmeantv  decay *nl  and  linear  operation  was  apptoa* he*l  \ 
value  **t  the  phase  shift  was  estimated  m this  opetatmg  iegu*n 
(**r  b*»(h  the  (>  ' an*l  I 0 11/  data  \t  frequencies  *»f  ' II  and 
above  the  wave  lot  ms  weie  sufficiently  sinusoidal  tv*  cxtta*l 
tellable  frequency  iesp**n\e  data 

I he  amplitude  latu's  and  phase  angles  **t  >0  >Al',  detei 
mined  fiom  the  time  histoty  data,  are  shown  plot  ed  as  a 


A 


function  i*l  excitation  trequcncv  far  tan  speeds  of  '560  and 
*-00  pun  m Fig  4 \l  each  speed  the  data  arv  compared  to 
I he  dequeue  v response  characteristics  ot  a simple  lirsl-ordet  lag 
having  the  torm 

&Q  hi  . 

AAI“,  itilfi.l 

I he  excellent  correlation  shown  on  (he  figure  indicates  (hat  a 
brst  order  function  adequatclv  describes  (he  dMtatnic  ctuiactei 
isdes  ol  (he  Ian  within  (he  (e>(  frequency  range  I tie  o*mer 
dequeue  v \x  appear*  (**  be  slightly  lower  at  (he  higher  Ian  spee*l 
1 he  value  ol  I , m I quadon  i I i i*  (he  inverse  slope  ot  (lie  tan 
pievsute  I low  chaiac  tertstk  at  (tie  mean  operating  point 


, I (ST  OAT  A 
. AH  A6K>  j Of 


»c  • 2 4 Hj 


FRIOUINCV  iM«» 

Fig  4a  Frequence  Koponsr  I Vila  it  2 <60  r pm 


FRIOOfNCV  (Mi* 

Fig  4b  h6|w>no  Response  hull  at  (.’TO  rpm 

Fig  4 - Bode  Plot  of  SFS  I00H  Scale  Model 
Fan  Response 

Dynamic  Analysis 

\n  analvsis  was  per(orme*l  t\*  identitv  the  principal 
dynamic  elements  of  (he  SI  S-IOOH  cushion  fan  model  s\  stent 
Included  in  (his  were  (he  following  tasks 

I Develop  an  analytical  model  which  includes  the 
principal  dynamic  elements  ol  the  tan  (es(  s\ stem 

Predict  (tie  behavior  which  each  ot  the  principal 
elements  would  produce  under  dvnamu  conditions 

( I xamine  and  compare  the  predicted  behaviot  with  that 
observed  during  the  test  program 

4 I Intimate  those  elements  which  do  not  produce  the 
observed  behavior 


< I xamine  in  greatci  detail  the  elements  whuh  tesult  in 
svstem  behavior  similar  to  that  observed  during  the  test  piogiam 
i|uan(if>  parameters  and  develop  picdklion  metlHvls 

\ schematic  diagiam  ot  the  tan  test  svstem  is  shown  in 
I ig  N I he  dvnamk  elements  which  are  shown  on  the  diagram 
are  il*  the  totational  inertia  and  torque  speed  properties  ol  the 
drive  motot  and  tan.  {2)  the  fluid  dviunik  properties  ot  tin  an 
within  the  tan  and  i *»  the  compliance  and  tl»*w  cturac tcnsticx 
ot  tlie  plenum  I hese  elements  weie  examined  as  possible  causes 
ot  the  observed  tan  behaviot 
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I ig  5 - Scale  M»h1oI  I an  Svstem  Schematic 

I he  equations  which  describe  the  svstem  aie  listed  below 
\s  listed  thev  are  in  a somewhat  smtplitied  torm  which  max  be 
expanded  in  more  detail,  it  icq  wired.  as  the  partieulai  elements 
are  examined 

System  Equations 

i See  I ig  < loi  Svstem  Schematic* 

Motor  Fan  Rotational  System 

I Motot  lorque 


1 1 S' 


Motoi  I an  Speed 


* I an  lorque 


I,  MS  OO 


Alternate  Solution 


I I an  Discharge  Pressure 


Alternate  Solution  twith  Ian  compliance  included! 

0,,  ft\  Al'( 


M,  .h0,,  -0, 


‘ "V  - 


2 Ean  Discharge  I Iom 

= I*.  H„ 

Plenum 

1 Mass  Accumulation 

\ *v«i  » 

2 Density 


3 Pressure 


4 Volume 


Z v, 


Aj  * effective  area  of  wall  section  i 
/,  wall  deflection  of  wall  section  i 


>V\  - m,  V»,V\A 


i i-’i 


1 1 )i 


I I4i 


,15l 


t I ol 


In  ill' termini-  t hi-  rtlivtj  ot  ipi'i-il  ilungn  I he  clla.li  ot  speed 
uiuliiMii  were  determined  Inun  Ihe  lolloitmg  hnean/ed  set  i*l 
dilution.  denied  from  I hr  list  nt  ei)rutiom 

J Motor  liirt|ue  (Irom  l'e|uatn>ii  t — 1| 

<tr„ 
di 


61. 


6 ij , 


r .’ill 


b I .m  l nri) ne  (Irnin  Iqualimi  iMi) 


ST. 


£i i 

llu'. 


6u.y 


31, 

iK7  *0, 


' l eakage  I lows 


c*p  - ( i>p  v 


: r. 


(17, 


list 


k,SceM  . k,6Q, 

i Motor- 1 .in  Rotation  (Irom  l i|iiation  , 5 r) 

'l^M  AI«  S'| 

d I an  Discharge  l*iessure  (neglecting  the  fan's  tin  id 
Jyiiaiiik  characteristics  1 1 v«m  Iquation  (M) 

3P,  <ll', 

M>-  ' 6u,“  * 5o7  SQ> 

R,  6uim  ♦ R, 

Substituting  l quations  \ 20)  and  t2lt  into  t22> 

Ks6Q.  0 


t:i » 


i 2 i 


1 .'41 


Ihe  tour  primary  differential  equations  which  describe  the 
dynamics  of  the  system  are  Equations  1 31,  1 1 1 ),  1 1 21.  and  1 13> 
Several  more  differential  equations  which  describe  the  wall 
motion  of  the  plenum  would  be  necessary  for  a complete 
description  However,  a knowledge  of  the  d\  namic  characteris- 
tics of  the  plenum  was  not  really  necessary  to  the  principal 
objective  of  the  study,  an  understanding  of  the  dynamic 
bchavior  ot  the  tan  Since  plenum  pressure  is  approximateh 
equal  to  the  tot al  pressure  rise  of  the  fan.  total  pressure  rise 
data  which  were  obtained  during  IOOH  testing  were  used 
directly  • l sing  the  pressure  rise  data  as  an  input  to  the  fan 
model  rather  than  a value  obtained  from  a simulated  plenum 
results  m a much  simpler,  more  accurate,  and  effective  was  of 
studying  fan  dynamic  characteristics 

Effects  of  Fan  Speed  Variations 

Ml  hough  speed  variations  w ere  not  regarded  as  a ser\  likely 
cause  of  tlie  fan's  dynamic  behavior,  there  were  some  questions 
as  to  the  abilits  of  the  speed  sensing  system  used  in  the  test  to 
record  dynamic  changes  therefore,  an  analysis  was  performed 


Ihe  following  matrix  equation  in  Laplace  notation,  formed 
from  t quations  t2.'t  and  i24).  describes  the  motor  fan  system 
with  tan  discharge  pressure  as  the  excitation  input 
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Solution  of  Equation  (25)  tor  fan  flow  iSQ,  in  terms  ot 
discharge  pressure  M’j  yields  ihe  following  transfer  function 
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Since  the  discharge  pressure  perturbation  5l*j  is  essentiallv 
fhe  same  as  the  variation  in  fan  piessure  rise,  the  latio  iSQ,  M’, 
is  the  same  function  evaluated  from  the  dy  namic  test  data  and 
plotted  as  a function  of  frequency  in  Lie  4 t herefore,  it 
variations  m fan  speed  were  the  cause  ot  the  tan  s dv  namic 
behavior  the  transfer  function  ot  Lquation  (2M.  evaluated  as  a 
function  ot  trequenev . should  vicld  a similar  function  An 
examination  ot  Equation  (2ot  and  its  parameters  indicates  that 


•It  should  be  pointed  out  (hat  the  purpose  ot  the  plenum  and  discharge  valve  was  to  provide  a variable  foukpteuutf  environment 
toi  the  tan  During  testing  Ian  flow  wav  observed  to  respond  in  a dvnamu  manner  to  the  oscillators  hackptessutc  As  the  equations 
indicate,  the  tan  will  respond  in  the  same  manner  to  a given  backprcsauie  vauation  regardless  of  how  the  variation  was  produced 
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wc  i or  net  frequency  trad  si 

I , l , 

lire  first  order  dvnamu  characteristic  ot  Fquation  t*M 
which  results  Irom  inertance  is  of  l he  form  indicated  h\  the  Ian 
test  data  us  plotted  in  I ig  4 In  I ouricr  lorm  Fquation  < >m 
becomes 


1 lt|-  I ♦ | 


in  I ig  4.  and  stall*  lan  dope  ol  - 2 ' It  s/lb  !t:  correspond  to 
an  inertance  value  0.0.1 1 5 tlb  it'  i 4 'though  t|u  reason 

lor  this  apparent  im tease  is  not  lulls  understood.  it  is  believed 
to  result  from  the  higher  flow  and  pressure  at  the  higher  speed 
I he  solute  is  more  unitormls  tilled  sstih  an  which  results  m an 
increase  in  effect ive  lloss  radius  and  path  length  I he  increased 
density  is  also  a contributor 

\n  evaluation  ot  the  lan  inertance  l , troin  l qualion  i I'l 
is  in  genoial.  sets  ditticull  1 he  tlow  path  is  recognized  as  verx 
wimples  and  the  lloss  itself  probjbls  not  uniform  and  unidtre*. 
tional  at  a given  o-ass  section  l herelore.  a reasonable  estimate 
ot  inertance  is  the  best  that  can  be  obtained 

It  the  tloss  ooss  section  \iSi  was  umionn  then  \ quation 
t *"*)  becomes  the  inertance  lor  a uniform  line  V units  m length 

■ PV  ’a  . , 
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Ohasht  performed  an  analytical  and  experimental  studs  i>t 
the  dsnanne  characteristics  ot  turbopumps  \n  analstical  model 
tor  the  frequenev  response  of  a fan.  derived  from  Ohashi’s 
model,  was  utilized  bv  lioldschmied  and  NNormles  4 In  the 
Ohashi  paper,  a model  tor  osullating  (low  through  a cascade  is 
developed  I’his  is  used  to  define  the  dvuamic  characteristics  of 
a turhopump  which  included  the  following  three  effects 

1 Oscillating  flow  through  a variable  area  conduit. 

2 Oscillating  flow  through  the  rotor  cascade,  and 
Oscillating  flow  through  the  stator  cascade 

Ihe  first  of  these  is  similar  tv'  the  inertia  effect  identified  as  the 
pnmars  contributor  to  the  fan's  d>  namic  properties  as  seen  m 
the  OINSRIK  test  ot  the  SI  s 100B  model  cushion  fan  Ihe 
latter  two  effects  are  mathematically  expressed  using  the  gamma 
function  and  the  C«auss  h\  pergeomelnc  function  Ohushi  goes 
on  tv>  present  a simplified  model  ot  a turbopump  which  considers 
only  the  rotor  and  stator  effects  I malls,  he  states  that  the  time 
constant  of  a turhopump  can  be  approximated  bs  a first-order 
lag.  which  is  the  same  tspe  ol  function  that  was  derived  bs 
analvsis  ot  the  IHVSRIM  test  data 

Ohashi’s  paper  concludes  with  a description  of  an  experi- 
ment to  measure  the  ds  namic  characteristics  of  a centrifugal 
pump  Ihe  comparison  between  his  simplified  theory  and  the 
experimental  data  showed  that  the  experimental  amplitude  ratio 
decreased  at  a faster  rate,  and  with  a larger  phase  shift  than  the 
theory  predicted  In  fact,  it  appears  that  a more  satisfactors  til 
between  the  experimental  data  and  theorx  could  be  made  with  a 
first-order  lag 

\ comparison  of  the  terms  ot  fquation  with  the 
numerical  values  determined  from  the  curve  fit  ol  the  test  data 
given  m Fig  4 requires  that 

-F,  *2.8  (40) 
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Now.  F,  was  defined  as  dQ,  . which  in  turn  was  evaluated 

from  the  25  M)  rpm  fan  characteristic  to  be  approximately 
-2  S fr  s/lb  ft*  Ihe  inertance  at  J200  rpm  appears  to  be 
somewhat  larger  \ corner  frequenev  F of  2 2 Hz.  as  indicated 


\ sketch  v't  the  tan  solute  is  shown  in  I ig  o Considering  the 
volute  only,  with  an  approximate  length  and  average  cross  section 
v'l  ^ 5 ft  and  I s ft*,  respectively  an  inertance  of 


(0  0024  KV  5) 
l 5 


0 0152  lb  s*  ft' 


is  oh  tamed  In  this  representation  the  tact  that  all  ot  the  flow 
does  not  pass  through  the  full  length  of  the  volute  tor  which  the 
inertance  sn.is  calculated  is  ottsol  bs  not  considering  the  inertance 
v'l  the  air  within  the  rotor  and  inlet 
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I ig  SI  S I00B  Scale  Model  Fan  l.oomeirs 

Xnv'ther  metluHl  of  estimating  the  fan  inertance  is  tv' 
determine  the  tloss  length  based  on  the  total  volume  of  the  fan 
and  an  average  flow  area  I he  total  volume  is  on  the  order  of 
20  ft’  l sing  an  average  flow  area  of  \ 5 ft-  the  inertance  i> 


ki  10  0024U20) 


0 021.1  lb  - v fts 


\ comparison  of  the  two  calculated  inertance  values  with 
the  value  indicated  bv  test  data  shows  them  tv'  be  ot  the  same 
general  magnitude  On  the  basis  of  this  relatively  close  agreement 
between  the  geometrically  calculated  values  and  those  determined 
Irom  test  data,  it  mas  be  concluded  that  air  inertia  is  a significant 
dvtiamu  element  and  is  the  principal  cause  of  the  observed 
dynamic  behavior 

I he  accurate  determination  of  the  inertance  of  a fan  which 
has  not  been  d\  namicallx  tested  is  difficult  as  prcvtouslx  mdiea 
ted  llv'wever.  for  prehminutv  svstem  analvsis  the  met  hints 


suggest  oil  will  give  an  adequate  representation  More  accurate 
values  can  then  he  obtained  through  dynamic  testing  as 
illustrated  here. 

Air  Compliance 

File  effects  of  the  compliance  of  the  air  within  the  fan  can 
he  evaluated  through  linearization  and  solution  of  fan  liquations 
(Si  through  1 12).  Since  the  length  of  the  How  path  through  the 
fan  (<  15  ft)  is  considerably  shorter  than  a quarter  wave  length 
(47.5  ft),  at  a frequency  of  6 11/  (the  highest  test  frequency)  a 
single  lumped  compliance  representation  is  adequate.  I he 
linearized  fan  characteristic  developed  previously  is 

6Q|:,  = (-  F;2)(6P,  - 6PJ  ) (45) 

This  relation  determines  the  How  into  the  fan.  I he  discharge 
How  from  the  fan  is  determined  by  the  inertance  of  the  volute 
as: 

6Q,  2 =<6P;.  - 6P,.)/Lf  (46) 

The  rate  of  change  of  mass  within  the  fan.  which  is  proportional 
to  the  difference  between  the  inlet  and  discharge  How.  is  then 
given  by 

6M,  = (60j , - 6Qj  ,)  (47) 

where  is  the  nominal  density  within  the  fan. 

Air  density  p,  within  the  fan  is  related  to  the  accumulated 
mass  of  air  M,  by 


and  similarly  for  small  changes: 

6p,  = 6M,  /V,  <4‘ 

where  the  volume  of  the  fan  Vj  . small  changes  in  fan  pressure 
6PJ  . and  density  6 p, , are  related  by 


6p;+p;  /sp.+^y 


Fig.  7 - Electrical  Analog  of  Lift  Fan  System 
including  Air  Compliance 


6Q,, 

P2 

M>, 

I.|  C's2  + (-  E , ) Lj  s + 1 

60, 2 

^T  = 

L,  C’s2  +(-F2)L,  s+  1 

For  the  SFS-100B  scale  model  fan.  the  volume  V,  is  20  ft3, 
and  C = 0.00644  ft5 /lb.  Using  values  of  F2  and  L(  previously 
determined,  the  functions  of  Equations  ( 53)  and  (54)  become 


(illi + ')  fe + ') 

fe + ')  (dr + ') 


and  in  linearized  form  by 

36  P!  7PI 

&K  = 6p.  = ^ Sp.  <51 » 

1 f*8pK  1 p,  1 

Taking  t he  derivatives  of  Equations  (49)  amt  (51 1 and  combining 
with  Equation  (47)  the  compliant  properties  of  the  fan  are 
described  as 

(SO  . -80,  , I 

6p;  = — r— 


Equations  (45).  (46).  and  (52)  then  describe  the  compliant 
and  inertial  characteristics  of  the  fan  to  small  perturbations 
These  three  equations  describe  a system  which  is  analogous  to 
the  simple  electrical  circuit  shown  in  Fig.  7 Solution  of  the 
linearized  equations  yields  the  following  transfer  functions  which 
give  the  How  rates  Q, , and  Q,  , as  a function  of  fan  discharge 
pressure. 


The  compliance  produces  the  high  frequency  form  at  421  rad  s 
in  Equations  (55)  ami  (56).  Ellis  frequency  is  much  higher  than 
the  lower  frequency  term  which  is  attributed  to  inertance  and. 
therefore,  is  of  negligible  significance  with  regard  to  the  fan's 
dynamic  characteristics  within  the  test  frequency  range. 

Digital  Computer  Solution  of  Fan  Performance 

A time  domain  simulation  was  developed  and  programmed 
for  solution  on  a digital  computer.  The  inputs  to  the  program 
include  the  steady  state  pressure-flow  characteristic  of  the  fan. 
the  fan  inertance  geometry,  and  the  discharge  pressure  time 
history  I'lie  simulation  utilized  a fourth-order  Kunge-Kutta 
integrator  to  integrate 


P]  = f(Q,  ) (determined  from  the  fan  characteristic) 
Pp  r P< I) (input  from  the  time  history  test  data) 


•Vx.m  * p;y 


! (the  calculated  density  of  the 
air  in  the  tan) 


Inertance  values  based  on  the  findings  ol  the  frequency 
response  evaluation  of  the  fan  and  the  dynamic  analysis  were 
used  in  the  simulation  Because  the  density  is  treated  as  a 
variable  within  the  nonlinear  simulation,  the  length  to  cross 
section  ratio  is  input  to  the  simulation  rather  than  the  inertance 
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value  as  indicated  in  I quations  (37)  and  (42).  A ratio  (V  A)  of 
75  ami  14  5 if  1 based  on  mcrtancc  values  of  0.0237  and 
0.0373  Ib-s*  ft*  were  simulated  for  the  fan  speeds  of  2560  and 
32(H)  i pm.  respectively 


Comparisons  of  simulated  response  and  test  data  at  a fan 
speed  of  2560  rpm  are  shown  in  I igs.  X through  13.  At 
frequencies  of  2 11/  and  above,  correlation  between  the  simulated 
and  recorded  How  response  is  excellent  Both  the  amplitude  of 
the  flow  perturbations  and  the  phase  angle  between  flow  and 
pressure  rise  are  in  close  agreement  At  these  frequencies  the 
perturbations  are  relatively  small  and  operation  remains  on  that 
portion  of  the  fan  characteristics  where  the  slope  is  nearly 
constant.  I he  svstem  behavior  is  normal  and  hence  the  simula- 
tion completed  describes  the  relationship  between  the  pressure 
and  flow  perturbations. 
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SIMULATION  Ml  SIJI  TS 


780  - 

1 740 
"r  .’oo  > 

S 160 

2 170 

80 


vf 

\! 


V*  // 

Y 

i — — i — 
0 v .10 
TIME  lMcon<U) 


\ - 


_i i I 


At  the  lowest  test  frequencies  (0.5  and  1.0  11/  shown  on 
Figs.  X and  0,  respectively)  the  rotation  of  the  discharge  butterfly 
valve  produces  very  large  pressure  and  flow  perturbations.  Under 
these  conditions  the  fan  is  operating  within  the  stall  region,  and 
at  cutoff  during  the  0.5  11/  perturbations,  for  a portion  of  the 
cycle  During  the  stall  portion  of  the  cycle  it  appears  that  the 
simulation  is  unable  to  determine  the  correct  pressure-flow 
relationship.  Although  the  flow  predictions  during  the  high  flow 
portion  of  the  cycle,  the  phase  relationships,  and  general  shape 
of  the  flow  history  appear  correct,  the  deep  troughs  within  the 


Fig.  10  - Correlation  of  Simulation  Response 
with  lest  Data  - Frequency  = 2 11/ 


H 


Fig  0 - Correlation  of  Siittulafiou  Response 
with  Test  Data  - Frequency  s I 11/ 


Fig  1 1 - Correlation  of  Simulation  Response 
with  lest  Data  - I- requeues  3 11/ 
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llow  trace  are  mu  achieved  lire  apparent  inadev|iiac>  inas  rest 
in  erllrei  lire  simulation  01  the  pressiue  use  data  which  weie  used 
as  the  simulation  excitation  Otiimg  the  tests  ol  the  SI  S lOOli 
svale  lair,  it  was  noted  that  the  rotation  ol  the  plenum  butletlls 


salve  was  not  steads  at  the  low  speeds  concspoiidiug  to  llow 
peituihation  lieipieimes  ol  U s ami  I 0 11/  heeause  at  these 
speeils  the  dine  ssstenr  developed  instil  In' rent  tonpie  to  lOimtei 
act  the  acuHlvitunin  immrents  which  were  acting  on  the  \al\e 


lig  I ' - Correlation  ol  Simulation  Response 
with  lest  Data  I requeues  s d 11/ 


l nder  these  conditions,  rapid  piessuie  llrn  tuatioiis  occ m r evl 
wuhin  the  plenum  and  at  the  tan  discharge  as  shown  in  I ig  S 
when*  a tlouhle  peaked  piessuie  transient  was  measuied  at  low 
llow s It  was  postulated  that  the  piessuie  tiaiisnut  Ural  occuned 
undei  these  conditions  was  actualls  nu»ie  seveie  than  was 
measuied  NVilli  stnh  an  exaggerated  piessuie  histoi\,  as  unluated 
m I ig  14.  the  simulated  llow  correlates  win  ilosels  with  test 
data 

It  is  also  possible  that  tan  ureitaine  \aires  as  a (unction  ol 
Ian  llow  rate  Such  a sanation  is  suggested  In  the  appaient 
d tide  re  nee  m \alues  observed  at  '>(»()  and  1 '00  rptu  luertaiice 
values  calculated  Hour  vis  nainw  test  data  obtained  at  the  two  tan 
speeils  ale  shown  plotted  as  a trmctioii  ol  llow  late  in  lag  I'' 

\ allies  deteimmed  Hour  estimates  ol  the  Ian’s  geomelis  aie  also 
show  n t oi  compaiisou  I xtiapolatiou  ol  the  data  obtained  at  tin- 
two  tan  speeds  would  indicate  that  vers  low  values  ol  ureitaiice 
mas  exist  at  low  llow  rates 

It  the  elleclise  meitance  vanes  directls  as  a function  ol 
llow,  then  the  evtiemelv  lapid  llow  changes  which  occui  m the 
stall  legion  would  be  possible  with  a less  prououneeil  pressure 
tlrntuatnm  than  indicated  in  I ig  Id  Mthough  a satiable 
meitance  was  not  simulated,  it  is  possible  that  such  a variation  is 
pa i (tails  lesponsible.  in  addition  to  the  pieviousls  discussed  eiroi 
m piessuie  use  data,  loi  the  disciepancs  in  low  tiequcncs 
simulated  (espouse  I he  variable  meitance  was  not  simulated 


I ig  I d C orrelation  ol  Simulation  Response 
with  lest  Data  I requeues  * t>  11/ 


I ig  14  Simulation  Response  with 
I sagger  a ted  Pressure  IVrlmbnlion 


NOMINAL  FLOW  (ft3/*) 

Fig  15  — Fan  Inerlancc  Variation  with  Air  Flow  Rate 

since  values  in  the  stall  region  would  have  been  difficult  to 
estimate  without  further  dynamic  testing  within  that  region. 

Also,  the  value  of  such  a simulation  with  regard  to  other  fans 
would  be  minimal  because  the  variation,  which  is  a function  of 
stall  and  flow  separation  characteristics,  would  undoubtedly  be 
different  for  every  far.,  if  it  does  exist. 

The  correlation  between  simulated  response  and  test  data 
indicates  that  the  single  lumped  inertance  representation  is 
adequate  for  lift  fan  dynamic  studies.  An  inertance  value 
estimated  from  fan  geometry  is  adequate  for  preliminary  studies. 
More  accurate  values,  possibly  as  a function  of  operating  condi- 
tions. could  be  obtained  from  dynamic  test  results. 

An  evaluation  similar  to  that  conducted  using  the  SFS-IOOB 
scale  cushion  fan  was  subsequently  conducted  at  DTNSRDC  on  a 
24-in.  diameter  centrifugal  fan  that  featured  a servo-controlled 
variable  geometry  <VC»)  device  to  control  the  fan  inlet  flow  and. 
ultimately,  the  plenum  pressure.  The  fan  was  a double-inlet, 
single  discharge  configuration  and  was  a 1 /4-scale  model  of  an 
early  design  of  the  lift  fan  for  the  2KSF.S.  Again  the  performance 
of  the  fan  as  measured  during  the  test  was  compared  with 
predictions  generated  using  a digital  computer  simulation  of  the 
tan.  In  this  case,  the  computer  simulation  was  expanded  so  that 
the  plenum  pressure  was  calculated  by  integrating  the  air  mass 


Fig.  16  - Correlation  of  Test  Data  and  Simulation  Response 
for  a Fixed  Geometry  Fan  Model  - Frequency  = 0.5  Hz 


accumulation  in  the  plenum,  rather  than  using  measured  test 
data.  A value  lor  the  fan  inertance  was  determined  for  the  fan 
volute  geometry  as  before,  and  in  addition,  inclusion  of  the 
inertance  of  the  plenum  discharge  valve  was  found  to  be  necessary 
to  achieve  proper  phase  correlation  between  the  valve  area 
oscillation  and  the  resultant  change  in  the  plenum  pressure. 
Results  of  the  comparison  of  plenum  pressure  test  data  and 
computer  predictions  are  shown  in  Figs.  16  through  19. 


Fig.  17  — Correlation  of  Test  Data  and  Simulation  Response 
for  a Fixed  Geometry  Fan  Model  - Frequency  = 4 H/ 


Fig.  IK  — Correlation  of  Test  Data  and  Simulation  Response 
for  a Fixed  Geometry  Fan  Model  - Frequency  = 0.5  Hz 


Fig.  19  — Correlation  of  Test  Data  and  Simulation  Response 
for  a Fixed  Geometry  Fan  Model  - Frequency  = 4 H/ 
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The  comparison  is  shown  at  How  oscillation  frequencies  of  0.5 
H/  and  4 Hr,  for  both  the  open  loop  (fixed  geometry  Ian)  and 
closed  loop  modes  of  fan  operation.  The  comparison  illustrates 
excellent  correlation  of  the  phase  response  during  open  loop 
operation.  The  simulation  slightly  overpredicted  the  peak 
pressure  at  low  frequencies;  at  higher  frequencies,  the  correlation 
of  amplitude  is  excellent.  When  the  fan  was  operated  in  the 
closed  loop  mode  with  variable  geometry,  the  same  comments 
apply  to  the  comparison.  In  these  cases,  the  simulation  included 
the  additional  state  variables  necessary  to  describe  the  motion  of 
the  VG  sleeve  devices. 

Conclusions 

The  conclusions  drawn  from  the  above  are  summarized  as 
follows; 

1.  Under  dynamic  air  How  conditions,  the  lift  fan  response 
8Q/5P  can  be  represented  by  the  transfer  function  of  a first- 
order  system  with  lag.  At  extremely  low  frequencies  on  the 
Bode  plot,  the  amplitude  ratio  is  flat  and  has  a value  correspond- 
ing to  the  negative  inverse  of  the  steady  state  fan  slope  dP/dQ. 
The  phase  angle  in  this  case  is  180  degrees.  As  the  frequency  is 
increased,  the  amplitude  ratio  decreases  and  the  phase  angle 
approaches  90  degrees.  The  corner  frequency  estimated  from 
the  test  data  was  somewhat  lower  at  the  higher  fan  speed  of  the 
test  data. 

2.  The  dominant  dynamic  property  responsible  for  the 
observed  behavior  of  the  fan  is  the  fluid  inertia,  or  inertance,  of 
the  air  within  the  fan,  A precise  evaluation  of  the  inertance  of 
a fan  from  purely  geometric  considerations  is  difficult  and 
analysis  of  test  data  is  required  for  an  accurate  determination  of 
the  inertance  properties.  For  most  dynamic  studies,  however, 
sufficiently  accurate  analyses  can  be  computed  from  fan 
geometry.  The  compliant  properties,  or  capacitance,  of  the  air 
within  the  test  fan  were  found  to  be  of  negligible  significance 
with  regard  to  the  observed  behavior  of  the  fan  under  dynamic 
flow  conditions.  This  is  also  believed  to  be  true  for  other  lift 
fans.  In  general,  these  effects  are  of  sufficient  magnitude  that 
they  should  be  considered  in  air  cushion  lift  system  dynamic 
analyses. 

3.  A digital  computer  program  has  been  developed  to 
predict  the  performance  of  lift  fans  operating  in  an  unsteady 
flow  environment.  Computer  predictions  of  the  fan  flow  rate 
have  been  correlated  with  measured  test  data  for  the  SFS-100 
scale  cushion  fan  and  a 1 /4-scale  model  of  a 2KSES  lift  fan. 
Excellent  agreement  between  the  predicted  and  measured  data 
is  evident  in  the  range  of  frequencies  investigated. 


The  work  reported  herein  was  performed  under  the 
sponsorship  of  the  Surface  Effect  Ship  Project  Office  (PMS-3041. 
Their  leadership  and  encouragement  are  gratefully  acknowledged. 
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